Quantum glass phases in the disordered Bose-Hubbard model 
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The phase diagram of the Bose-Hubbard model in the presence of off-diagonal disorder is deter- 
mined using Quantum Monte Carlo simulations. A sequence of quantum glass phases intervene at 
the interface between the Mott insulating and the Superfluid phases of the clean system. In addition 
to the standard Bose glass phase, the coexistence of gapless and gapped regions close to the Mott 
insulating phase leads to a novel Mott glass regime which is incompressible yet gapless. Numerical 
evidence for the properties of these phases is given in terms of global (compressibility, superfluid 
stiffness) and local (compressibility, momentum distribution) observables. 
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The competition between disorder, interactions and 
commensurability in quantum many-body systems is 
known to produce novel quantum glassy phases, char- 
acterized by a gapless spectrum and by the absence 
of a global order parameter [ij. ''He adsorbed on 
porous media0|, granular superconductors [1], disordered 
magnets Q are but a few manifestations of localization ef- 
fects due to random potentials in interacting bosons. In 
recent years, ultracold atomic gases in magneto-optical 
traps have opened a new frontier in the study of strongly 
correlated systems, as unprecedented control over exper- 
imental parameters in these systems makes them ide- 
ally suited for studying many-body phenomena. Dis- 
order can be generated in optical lattices by exposure 
to speckle lasers 0, incommensurate lattice-forming 
lasers 0, H, 0] , and by other means [l^ . The interplay be- 
tween disorder and interactions in trapped Bose-Einstcin 
condensates has recently been explored experimentally in 
^^Rb, both in the continuumj^, [ll|, [13] and in an optical 
latticed, [H. 

Theoretically, the effects of random potentials on inter- 
acting bosons in p eriodic lattices have been studied using 
analytic^, [3, E, [3] and numerical techniques [l3|. It is 
now well established that even infinitcsimally small po- 
tential disorder can destroy the direct superfluid (SF) 
to Mott insulator (MI) transition in one dimension by 
introducing an intervening insulating, but compressible, 
Bose glass (BG) phase[ll, [3]. Surprisingly, while the 
effects of potential disorder have been widely investi- 
gated, other kinds of disorder, e.g. hopping or inter- 
action strengths, have remained largely unexplored until 
recently, when it was demonstrated that these can lead 
to qualitatively quite different phenomena [l^, 0, 21 1. 



In the quantum rotor model with off-diagonal disorder, 
the SF to MI transition takes place via an intermediate 
Mott glass (MG) phase 0,[2l|- a unique incompressible, 
yet gapless, glassy regime that was first reported in disor- 
dered fermions with extended range interactions [l^. It is 
thus of great interest to investigate whether such a phase 
also appears in the Bose-Hubbard model (BHM), and if 
particle-hole symmetry is essential for the stabilization 
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FIG. 1: (Color online) Phase diagram of the 1-D BHM with 
spatially correlated off-diagonal disorder: = t/3 and t 

with equal probability. For comparison, the phase boundaries 
of the clean system (fi,i+i = t) are shown with dashed lines. 
The phases are characterized by the stiffness ps and global 
compressibility k: superfluid (SF) [ps > 0, k > 0), Bose glass 
(BG) {ps = 0,K > 0), Mott insulator (MI) {ps = 0, k = 0). 
The Mott glass (MG) and MI phases have identical global 
properties and differ only in local properties. 



of such a MG. 

We use large-scale QMC simulations to study the ef- 
fects of off-diagonal disorder in the Bose-Hubbard model 
on a one-dimensional (ID) lattice. We find that, in con- 
trast to diagonal disorder, the Mott lobes do not shrink 
in the presence of off-diagonal disorder. Instead, there is 
an extended BG phase separating the MI lobes from the 
SF regime, and an additional MG regime emerges. 

The Bose-Hubbard model is given by the Hamiltonian 



u 



(1) 

where &| ( bi) creates (annihilates) a boson at site i, i^.i+i 



is the hopping matrix element between sites i and i + 1, 
U is the on-site interaction strength, and /x is the chem- 
ical potential. OfF-diagonal disorder is introduced in the 
form of a bimodal distribution of the hopping matrix ele- 
ments - = t or t/3 with equal probability 23[. The 
Stochastic Series Expansion (SSE)[2J| method is used to 
simulate the model([T]) in chains of length L = 64 — 256. 
For each set of parameters, 200-1000 realizations of dis- 
order are sampled. To characterize the emerging phases, 
we compute the superfluid stiffness and the global com- 
pressibility. In simulations employing updates that sam- 
ple different winding number sectors, the stiffness is con- 
veniently obtained from the fluctuations in the winding 
numbers of the world lines as Ps = [{W'^) /2f3]av, where 
[■■■]av denotes averaging over multiple realizations of dis- 
order. The global compressibility, k, is the energy cost 
of adding a particle to the system. It is defined by 
K = — {n)'^]av, where n is the density of particles. 

At sufficiently large U, local variations of the hop- 
ping amplitudes create a non-trivial landscape of higher- 
mobility domains ( "lakes" ) that coexist with gapped, lo- 
calized regions. In order to tunc the size of these lakes, 
we introduce spatially correlated disorder. The Fourier 
filtering method [IHl is used to generate a sequence of ran- 
dom numbers, r](i), with long-range algebraic correlation. 



C(z-j)^Mz),7(j))~|z-j|-". 



(2) 



Without any loss of generality, we choose a = 0.3. This is 
then mapped onto a bimodal distribution of hopping inte- 
grals, ti^i^i, between sites i and i + 1 of a one-dimensional 
lattice. While the qualitative features are same for cor- 
related and random disorder, spatially correlated disor- 
der favors extended domains in the disorder realizations. 
This, in turn, allows for a more reliable identification of 
the phase boundaries of the unique MG phase which dif- 
fers from the MI only in terms of short-range properties. 
Its proper characterization thus depends crucially on the 
sampling of finite-ranged domains. With uncorrelated 
disorder, typical domain sizes are smaller, which renders 
the MG phase more difficult to detect numerically. 

The results of the simulations are summarized in Fig.[Tl 
rendering a rich phase diagram. In the absence of disor- 
der, the ground state of the BHM is in either the SF or 
the MI phase. The SF-MI transitions (dashed lines) are 
marked by the simultaneous vanishing of the compress- 
ibility, K, and stiffness, ps, leading to the well-known se- 
quence of MI lobes, with integer fillings 2^, 27|, as p/U is 
varied. In the presence of diagonal disorder, the two tran- 
sitions - Ps = and K = - decouple and a compressible 
(k > 0), insulating {ps = 0) BG phase appears. With off- 
diagonal disorder, an additional glassy phase ~ the MG - 
is realized. The MG has global properties identical to the 
MI {ps = 0, K = and integer filling), but differs in local 
properties. In contrast to diagonal disorder, the Mott 
lobes do not shrink with respect to the clean case. There 
are no glassy phases in the atomic limit {t/U — > 0). At 




FIG. 2: (Color online) Stiffness (ps) and compressibility (k) as 
functions of the chemical potential (p) at constant t/U =0.16 
(cut along the dashed line in Fig[T|). For comparison, data of 
the clean system are also shown. In the absence of disorder, 
stiffness and compressibility vanish simultaneously at the SF- 
MI boundary. In the presence of disorder, the pa — and the 
n — transitions decouple, resulting in the glassy phases BG 
(ps = 0, K > 0) and MG {ps = 0, k = 0). The MG regimes I, 
II and III differ from the MI in terms of local quantities. Note 
that there is a finite MG phase with n — 2 (1.28 < n/U < 
1.34) while there is no corresponding MI phase. 



finite t/U, the direct SF-MI transitions of the clean limit 
are replaced by SF-BG-MG-MI sequences. Our results 
are consistent with all the transitions being continuous. 

Details of the simulations are illustrated in Fig. [21 
The global stiffness and compressibility are shown as 
a function of the chemical potential p/U at constant 
U/t ~ 6.25, i.e. along the dashed line in Fig.[TJ The data 
for the clean system is shown for comparison. At small 
values of p/U, the ground state is SF with ps > 0, k > 0. 
With increasing p/U, the system passes through a series 
of phases, as indicated in the figure. The first SF-BG 
transition is marked by the vanishing of ps, while k re- 
mains finite across the transition. This coincides with the 
SF-MI boundary in the clean system with ti^i+i — t/3 for 
all links. In the BG phase, the ground state is a mixture 
of SF and MI domains of all sizes - the domains with 
ti,i+i = t/3 are local Mis, while those with i^.i+i = t 
are locally SF. The MI phase of the disordered system 
coincides with that in the clean limit with i^^i+i = t. Of 
particular interest are the regions marked I, II and III. 
For these ranges of p/U, the global compressibility, k, 
and the stiffness, ps, vanish identically. The ground state 
has integer filling (n = 1 for regions I and II and n = 2 
for region III). The global properties are thus identical 
to those of the MI phase. As shown next, unlike the MI 
phase, the ground state has locally compressible regions, 
and these ranges of p/U can accordingly be identified 
as MG phases. Since these regimes lie outside the Mott 



lobes of the clean system with ti^i+i = t (in particular, 
there is no n = 2 MI phase for this value of U/t), the 
ground state consists of a mixture of SF = t) and 

MI (^i.i+i = i/3) domains as in the BG phase, but differs 
from the BG by being globally incompressible. 

Having explored the global properties, we focus on the 
local compressibility and the momentum profile of the 
ground state in each phase. The local compressibility 
at site i is defined as the local number fluctuation, Kj = 
/3[(n|) — {ni)''^]av, where Ui is the particle density at site i. 
The momentum distribution is obtained from the equal- 
time Green's function. 
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n{q 2ti) measures the short-range coherence in the 
system. The left panel of Fig. [3] shows the distribution of 
the local compressibilities in a lattice of length = 128, 
averaged over 800 disorder realizations. In the MI phase, 
the distribution is peaked at small Ki . Conversely, in the 
SF phase it is peaked at a finite value. In the BG and 
MG phases, the distribution has a double-peaked struc- 
ture, consistent with co-existing SF and MI domains in 
these phases [2^. Fig. [3]Jb), shows momentum profiles in 
these phases. n{q) is sharply peaked at g = in the SF 
phase and exhibits only a weak broad maximum around 
(7 = in the MI phase. The behavior of the MG and 
BG phases is similar and intermediate between the two 
limits - there is a peak at g = arising from the SF 
domains, but the height is reduced due to the MI do- 
mains. Thus the BG and MG phases have very similar 
short-range structure (both phases consist of co-existing 
SF and MI domains), although their global responses are 
rather different, as documented by Fig. [51 While the 
MG phase is globally incompressible, there exist a gap- 
less channel that allows the exchange of particles between 
adjacent SF and MI domains, leading to its incompress- 
ible yet gapless character. Such a mode is made possible 
due to the local variation of the kinetic energy of the 
bosons in the different domains. This explains why the 
MG phase is not observed for purely diagonal disorder 
and is unique to off-diagonal disorder [2^. The need to 
probe local properties also makes it difficult to detect this 
phase in current optical lattice experiments, where, ad- 
ditionally, domains of different phases co-exist. Current 
developments in local probes using selective microwave 
spectroscopyjs^ and the recently suggested method of 
using induced controlled interactions between a "probe 
particle" and the many-body state appear promising (3ll|. 

To further characterize the various phases and probe 
the gapless nature of the MG phase, we study the static 
structure factor, S{q), for the density-density correlation 



S{q) 



■rk), 



(4) 




1000 



100 



0.01 0.1 
q/27t 



FIG. 3: (Color online) (a) The distribution of local (site) com- 
pressibilities in the different phases averaged over 800 realiza- 
tions of disorder. The double-peaked distribution in the BG 
and MG phases indicate the presence of co-existing SF and 
MI domains, (b) Momentum profiles of the corresponding 
ground states. The BG and MG phases have similar momen- 
tum profiles, reflecting the similarity in their local structure. 



and the associated susceptibility. 
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A finite value of S{q) as g ^ implies the phase is com- 
pressible, i.e., the energy gap towards adding or removing 
a particle vanishes in the thermodynamic limit. Con- 
versely, a vanishing S{q) as g — > is a signature of an 
incompressible phase with a gapped charge excitations. 
Additionally, the ratio 2x{q)/S{q) as (7 — > gives an up- 
per bound for the charge excitation gap (3^. Fig. [4] shows 
the structure factor and the susceptibility as a function 
of momentum for representative points in each of the 
different phases. Consistent with global compressibility 
measurements, S{q) vanishes as q — > in the MI and MG 
phases, whereas it remains finite for SF and BG ground 
states. The upper bound of the charge excitation gap, 
A = 2x{q)/S{q) as q ^ 0, on the other hand, sheds light 
on the nature of the MG phase. Although Amg is finite, 
it is less than the effective chemical potential difference 
^f^Aic/ between type I (t = 1) and type II {t = 1/3) do- 
mains. While there is a gap towards overall addition or 
removal of bosons, there exists a gapless mode involving 
the transfer of particle across the boundaries of the two 
types of domains which is driven by the effective potential 
energy difference between adjacent sites at the domain 
boundaries. This is in contrast to the MI phase where 
the effective potential difference S^^^j^ between the two 
types of domains is less than Ami- This confirms the 
existence of a gapless mode in the MG phase. 

In conclusion, using large-scale quantum Monte Carlo 



4 




I I 


(b) 


. „ I II 






iT e III 
' OLl - 

/ MI 


1,1,1, 



05ER46240 (SH). S.H. ackowledgcs the hospitality of the 
Los Alamos National Laboratory where part of this work 
was carried out. The simulations were carried out at the 
high-performance computing computing center at USC. 



X 







0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 
q/27i q/27t 

FIG. 4: (Color online) (a) Static structure factor, and (b) 
static susceptibility for the diagonal correlation in the differ- 
ent phases as a function of momentum. The effective poten- 
tial difference between type \ (t = V) and type II (t = f/3) 
domains in the MG and MI phases are denoted by 
and respectively. Consistent with global compress- 

ibility measurements, S{q) vanishes as q — > in the MI and 
MG phases, but remains finite in the SF and BG phases. 
Aj\/G > 2x(g)/'S'(g) as g — > in the MG phase, implying the 
existence of a gapless mode. 



simulations, we have observed that off-diagonal disorder 
(random hopping), yields a non-trivial sequence of quan- 
tum glass phases that intervene at the interface between 
the Mott insulating and the superfluid regimes of the 
clean system. In particular, the coexistence of gapless 
and gapped regions close to the Mott insulating phase 
leads to a novel Mott glass phase which is incompressible 
and gapless. It shares some of the global properties of 
the Mott insulator, but resemble the Bose glass in local 
properties. It is remarkable that the phase boundaries 
of the Mott insulating phase are basically unaffected by 
ofF-diagonal disorder, and that the extent of the Mott 
glass can be tuned by varying the range of the spatial 
correlations in the disorder realization. It is also evident 
that particle-hole symmetry is not a pre-requisite for ob- 
serving the Mott glass phase. The essential ingredients 
for its realization are off-diagonal (hopping) disorder and 
commensurate filling. Finally we observe that while the 
present results are obtained for a one-dimensional sys- 
tem, the Mott glass phase is expected to be more ro- 
bust in higher dimensions. As seen here, the robust- 
ness of the Mott glass phase depends crucially on hav- 
ing large domains with uniform intra-domain hopping. 
This is achieved more readily in higher dimensions due 
to increased co-ordination number. Indeed, in three di- 
mensions, the Mott glass phase is expected to persist to 
finite temperatures. 
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